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Abstract 

We consider the asymptotic behavior of the incomplete gamma functions 7( — ffl, —z) and 
r(— a, —z) as a — > oo. Uniform expansions are needed to describe the transition area z ~ a, 
in which case error functions are used as main approximants. We use integral representations 
of the incomplete gamma functions and derive a uniform expansion by applying techniques 
used for the existing uniform expansions for 7(ffl, z) and T(a, z). The result is compared with 
Olver's uniform expansion for the generalized exponential integral. A numerical verification of 
the expansion is given in a final section. 
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1. Introduction 

The incomplete gamma functions are denned by the integrals 

t*z poo 

j(a,z)= t^e-'dt, T(a,z)= f^e^dt, (1.1) 

JO J z 

where a and z are complex parameters and t a takes its principal value. For 7(0, z) we need the 
condition > 0, for T(a,z) we assume that | argz| < n. Analytic continuation can be based 
on these integrals, or on series representations of 7(0, z). We have 7(0, z) + T(a, z) = T(a). 
Another important function is defined by 

7>,z) = ^ 7 («^) = ^T fu^e—du. (1.2) 

This function is a single-valued entire function of both a and z, and is real for positive and 
negative values of a and z. For T(a, z) we have the additional integral representation 

e~ z f°° e~ zt t~ a 
T(a,z) = — i r/ 7— dt, »a<l, &z>0, (1.3) 
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which can be verified by differentiating the right-hand side with respect to z. 

For other information on the incomplete gamma functions we refer to Chapter IX of the 
Bateman Manuscript Project (1953) or Chapter 11 of Temme (1996). 

The purpose of the paper is to derive new uniform asymptotic expansions for the functions 
7(— a, — z), T(— a, — z) and7*(— a, — z). The uniform expansions for 7(0,2) and T(a,z) of our 
earlier papers, which will be summarized in the next section, are not valid for negative values 
of a. 

Recently there is much interest in asymptotic properties of incomplete gamma functions. 
The function T(a,z), or the related exponential integral 



plays an important role in Berry's smooth interpretation of the Stokes phenomena for certain 
integrals and special functions; see Berry (1989). Olver (1991a) investigated E p (z) in 
particular at the Stokes lines argz = ±7r and used the results in Olver (1991b); see also 
Olver (1994). We summarize Olver's result in the next section. In Dunster (1994a) a new 
expansion for the function E p (z) is given for complex values of z and large positive values of 
p and in Dunster (1994b) error bounds are given, in particular for the Stokes' smoothing 
approximations. In Paris (1994) the uniform asymptotic expansions of the next section are 
used for complex values of the parameters in a new algorithm for computing the Riemann zeta 
function on the critical line. 

The computational problem of the incomplete gamma functions for complex values of a 
and z is not well solved in the software literature. In particular when the complex parameters 
have large negative real parts existing computer programs may give really false answers. We 
expect that the new expansions for 7(— a, — z), T(— a, — z), 7* (— a, — z) and Olver's expansion 
for Ep(—z) will be of value to solve this problem. 

2. Uniform expansions for incomplete gamma functions. 

We summarize the known uniform expansions for the incomplete gamma functions; see Temme 
(1975), (1979) and more recently Olver (1991a), (1994). 

2.1. Uniform expansions for P and Q. 
Let r] be the real number defined by 



Extend the relation between A and rj to complex values by analytic continuation. We have 




(1.4) 



\r\ 2 = A — 1 — In A, A > 0, sign(??) = sign(A — 1). 



(2.1) 



J ?~(A-l)[l-i(A-l) + ...] J A-l. 



Then we write, with A = z/a, 




(2.2) 



P{a,z) 



T(a) 



\ eik(-riy/a~/2) - R a (r]). 
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The error functions are denned by 

e _t dt (2.3) 



erf z = [ e 1 dt, erfc z = 1 — erf z = —=■ [ 
v 71 " Jo v 71 " ■/* 



The error functions are the dominant terms in (2.2) as a tends to infinity, and they describe 
the transition at a = z. 
We have 

R a ( v ) = -==S a ( V ), 5 a (r?)~^^, (2.4) 
-vZttci — 1 ft 

as a — > oo, uniformly with respect to | arga| < tt — 6\, and A, | argA| < 2n — 62, where 61,62 
are arbitrarily small positive constants. More information on the coefficients C n {rj) is given in 
the next subsection. 

In this paper we derive similar new expansions for 7 (a, z) and T(a, z), where we concentrate 
on negative values of a and z. Because simple reflection formulas are missing that relate 7(— a, z) 
with ■y(a,z), etc., we cannot simply transform the above expansions to this case. 

2.2. Further details on the previous results. 

Because the new expansions (see §3) are closely related with the one in (2.4), we give some 
details on computing the coefficients and on the mapping A — > 77(A) defined in (2.1). 

The expansion in (2.4) for S a (rj) can be obtained by differentiating one of the equations 
in (2.2) with respect to 77, which gives 

±S a ( V ) - a V S a ( V ) = a [1 - f(v)/T*(a)] , 

where 

1 dX 



HV) = = r» = <yA7(2^) e a a- a T(a). (2.5) 

Substituting the asymptotic expansion for S a (r]) one finds for the coefficients the relations 

0,(77) = —^---, V C n (v) = C n _ 1 (77)+ 7n /(r7), n>l, (2.6) 
A — 1 77 drj 

where 7 n are the coefficients in the reciprocal gamma function expansion 

1 In 



T*(a) ^ a 

V ' n = 



a ^°°; |arga|<vr, (2.7) 



of which the first few are given by 



139 571 



70 = 1, 71 = - — , 72 = — , 73 = eiOAn , 74 



12' ' 288' ' 51840' ' 2488320 

Next we summarize the properties of the mapping A — > 77(A) defined in (2.1). More details 
are given in Temme (1979). Write 77 = a + i/3. The mapping is one-to-one for — 2n < arg A < 
2ir, A 7^ 0. The corresponding 77— domain is given in Figure 1. The mapping is singular at 
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Figure 1. Corresponding points in the A— and ^—planes for the mapping A — > r](X) defined 
in (2.1). The points B ± are singular points in the rj— plane. In the A— plane they correspond to 
exp(±27ri). The ^—points D ± 7 E ± correspond to A— points D ± ,E ± having phases ±ni. The 
thick branches of the hyperbolas are branch cuts. The parabola shaped curve in the A— plane 
corresponds to the imaginary axis in the rj— plane. 
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A = exp(±27ri) with corresponding points B ± = ^/2k (—1 ±i). The positive real A— axis (with 
argA = ±27r) is mapped to branches of the hyperbolas a(3 = ±27r, a < —\p2/K. When we 
cut the 77— plane along the two thick parts of the hyperbolas in Figure 1 we have a one-to-one 
interpretation of the mapping in (2.1), and the function 77(A) is singular at 77 = y/2~K (—1 ± i). 
The parabola shaped curve in the A— plane is the set of points 9(A — In A — 1) = 0, and given 
by the equation p{0) = 0/ sm6, —tt < < tt, p(0) = 1, where p, are the polar coordinates of 
A : A = pexp(i#). The corresponding set in the 77— plane is the imaginary axis. 

2.3. Olver's expansion. 

Olver (1991a) investigated the generalized exponential integral E p (z) defined in (1.4). Let 

py 1 2vr zp- 1 
and z = pe l6 , a = n — p + p with p a large parameter, p fixed. Then 

■\p{c{6)Y 00 



F n+p (z) ~ {-l) n ie 



ierfc \c.(0)^\p ] - l E ^_ Y, (s). «) ( J ) 



(2.8) 



uniformly with respect to 6 G [— n + 3n — 6] and bounded values of \a\; 6 denotes an arbitrarily 
small positive constant. Furthermore, 



c(0) = ^2 {e ie + i(0 - tt) + 1} , 

with the choice of branch of the square root that implies c(0) ~ (n — 0) as ^ — > 7r; the coefficients 
52s(^,«) are continuous functions of ^ and a. A similar expansion for F n+P (z) is given when 
G [— 37T + 7r — £]. As Olver remarks, this expansion quantifies the Stokes phenomenon, that 
is, the rapid but smooth change in form of other expansions as passes through the common 
interval of validity of the other expansions. 

In the present paper we give an expansion that is also valid in a domain around the 
negative z— axis and that also contains the error function. Our expansion is quite different, 
however, and is strongly related to the uniform expansion given in (2.1)-(2.4). 

3. Starting from negative values. 

We start with (1.3) and replace a with —a: 

-z poo -Zt iO 

T(-a, z) = — — / — — - dt, 5Ra > -1, Mz > 0. (3.1) 
T(l + a)J t + l 

By turning the path of integration and invoking the principle of analytic continuation we can 
enlarge the domain of z. For example, when we consider the path in (3.1) along the positive 
imaginary axis, the integral is defined for — tt < arg z < 0, and in the overlapping domain 
— 7r/2 < arg z < its value is the same as in (3.1). We turn the path from the positive 
imaginary axis argt = ir/2 to the negative axis argt = tt, avoiding the pole at t = —1 by 
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using a small semi-circle. The integral is then defined for — |7r < argz < — ^tt. We change the 
variable of integration and the result is: 

e Zgi7ra poo g-zt^a ^ \ 

T(-a,ze~ l7T ) = — -/ dt, > -1, --ir < arg z < -it (3.2) 

V ' ; T(l + a) J t - 1 ' ' 2 B 2 v ; 

and we avoid the pole by integrating under the pole at t = 1. By turning the path in (3.1) 
clockwise we obtain: 

e z e~ iva f°° e~ zt t a 1 1 
T(-a,ze +iv ) = — -/ dt, »a>-l. 7r<arg^<-7r (3.3) 

V ' ; r(l + a) J t - 1 ' 2 8 2 V ; 

and we avoid the pole by integrating above the pole at t = 1. An easy consequence is 

_27T? 

r(-a, ze + ™) - e-™ a T(-a, ze~™) = - (3.4) 

T(a + 1) 

which follows from computing the residue of the integral over the full circle around the pole at 
t = 1. It also follows from the fact that T(a,z) = T(a)[l — z a ^*{a,z)] and ^*{a,z) is entire. 

We proceed with (3.3) and assume, for the time being, that a and z are positive. We 
change the variable of integration t — > at/ 'z, and obtain 

zx - a iav f00 dt 



z\-a -ta-K poo 

T (-a, ze™) = , / e~ at t c 

{ ' r(« + i) J 



t-x 



„\an 2 —ia-K roc i/- ' 

= TFTTT / ^VC)^ 
r(a + i) c-77 

where r\ = tj(X) is defined in (2.1), A = z/a, ( = r](t), that is, 

=t-\nt-l, signC = sign(i - 1), 

and 

5(c) = dc^A = r^i^A- (3 - 5) 

In the ^— integral the path passes above the pole at ( = rj. 

When z ~ a, that is, A ~ 1, the pole is near the saddle point ( = 0, and for large values 
of a we need an error function to describe the asymptotic behavior. We can split off the 
pole by writing g(Q = [g(Q — g(v)] + 9(v), where, as is easily verified, g(rj) = 1. Using the 
representation 

e z ' f°° e"* 2 
^eiic iz = / dt z 6 C, 

2 2™ y.^ t - Z 

where the path passes above the pole at t = z, we obtain (introducing a suitable normalization) 

-r(a + 1) — r (-a, ze**) = ±erfc Wf - ,-—^(7?), (3.6) 
^7r« v v / V27ra 
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and, for instance by using the relation in (3.4), 

r(a + 1} r (_a ' ze ~^ = ^ erfc i~ iv \f%) + i S= Ta(??) ' (3 - 7) 



where 



T.(,) = - 1 /ijr«-*< , »K)«, W<) = ^W. (3.8) 



Standard methods for integrals can be used now (see Wong (1989)) to obtain an asymptotic 
expansion of T a (r]) in negative powers of a. It is easier, however, to use a differential equation 
satisfied by T a (rj). In this way we can identify the coefficients of the expansion with those in 
(2.4). 

Differentiating (3.6) with respect to 77, and using (1.1), (2.1) and (2.3) we obtain the 
differential equation 

±T a ( V ) + ar,T a (r,) = a [f( V )T*(a) - 1] , (3.9) 

where f(rj) and T*(a) are given in (2.5). As in (2.7) we have (with the same coefficients 7„ 
indeed) 

oo 

r *( ffl )~£(-l)"|l, a^oo, |arga|<7r. 

n = 

Substituting this and the expansion 

t („) ~ f>ir^M (3.io) 

n = 

in (3.9) we find for C n (rj) the same recursion as in (2.6), and we conclude that in (2.4) and 
(3.10) the same coefficients occur. 

3.1. The domain of validity of expansion (3.10). 

We return to the case that a, z and A = z/a are complex parameters and claim that expansion 
(3.10) holds as a — > oo, uniformly with respect to | arg a\ < tt — 6± and | arg A| < 2ir — 62, with 
^1,^2 arbitrarily small positive constants. 

To verify this we observe that 77(A) is analytic and univalent in the sector | arg A| < 2n. 
This follows from §2.2. The function ((t) is in fact 7](t) and the quantity t = t(Q occurring in 
g in (3.5) is singular at the points B ± shown in the lower part of Figure 1. The functions g 
and h (given in (3.8)) are singular in the same points. It follows that h(Q is analytic in the 
disk I C| < 2y/n and in the sectors < When a is complex we can turn the path of 

integration in (3.8) inside the sectors \^s(\ < and 3?a( 2 can be kept positive on the path 
of integration as long as | arga| < tt. This gives the domains for a and A. 

We remark that the function h(() is bounded at infinity inside the sectors < |3?C|- 
This property is not needed to verify the above result, however. 
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3.2. A real expansion for 7*( — a, —z). 



It is of interest to have a result for the function 7* (a, 2;) (see (1.2)) for negative values of the 
parameters. In that case 7* (a, z) is real and we have a real asymptotic representation: 



7*(-a, -2) = z 



cos 7ra 



2a e |£«J 2 



sin7ra 



i^(Vt) + ^ (r?) ]}' (3 - n) 



where 



/" 

JO 



e* di 



which is Dawson's integral. The relation with the error function is: 

F(z) = — ji^n e~ z erfiz. 
Representation (3.11) can be verified by using (3.3) and the relations 

T(-a, ze +iv )T{-a) [l - ^e^V (-a, -z)} , T(a + l)T(-a) 



IT 



sin7ra 



In (3.11) we see that the oscillatory behavior of 7*(— a, —z) is described by two terms, one 
with cos7ra and one with exp(-^arj 2 ) sin7ra (the factor containing Dawson's integral is slowly 
varying when the parameters a and z are positive). This complicated oscillatory behavior 
is one of the problems in writing reliable software for the functions 7(11, z) and T(a, z) when 
the parameters have large negative real parts. Dawson's integral becomes dominant when we 
consider complex values of the parameters. The function F(z) is an entire odd function and 
has for $lz > the asymptotic behavior 



F(z) ~ < 



2z 



if argz < \ir, 



00, 



^ sign(^2;)i«V7r e~ z elsewhere. 



3.3. Comparing Olver's expansion with our expansion. 

In Olver's expansion (2.9) the modulus of z is the large parameter, whereas in our expansion 
(3.10) a is the large parameter. In addition, in Olver's expansion the transition occurs at 6 = 7r, 
whereas in our expansion the transition takes place at 7] = 0, that is, at A = 1, with A = z/a. 
We have r\ ~ A — 1, as A — > 1. Hence, when z = a + ie, then r/ ~ ie/a. When a is positive 
and e increases from negative to positive values, ze* % crosses the negative axis with increasing 
argument, and ^erfc(«ry^/ a/2 ) of (3.6) changes rapidly from to 1. 

Olver's expansion has more parameters than ours because his analysis started by consider- 
ing F n+P (z) as a remainder in the expansion for F p (z). However, one can take n = 0,p = a + 1. 
Then Olver's notation agrees better with our notation. In his result the parameter a = p — p 
is assumed to be bounded, and a restriction like this is not present in our method. 
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Remark. We verify what happens when the parameters a, z in (2.2) are taken with negative 
signs. First let 3a, ^sz be positive, and replace in Q(a, z) the parameters a, z with ae _7ri , ze -7r % 
respectively. The quantity 7] defined in (2.1) does not change by this operation, whereas the 
expansion for S a (r]) becomes the expansion for T a {rj) given in (3.10). When we formally write 
S- a (rj) = T a (r]), which certainly is not true, it follows that the right-hand side of the first line 
in (2.2) formally becomes 



ierfc -ivVa/2 + i-=T a (ri) 



which is the right-hand side of (3.7). But 

T (—a ze~ 7ri ) 1 

Q(-a, z~* 1 ) = v / — - — '- = T(a + 1) sinvrar(-a, ze~ vt ), 

T( — a) tt 

which means that the left-hand side of (3.7) equals 

T(a + 1) T (-a, ze-™) = \— Q (-a, ze'™) . 

Hence, by proceeding formally from the relation for Q in formula (2.2), we miss the factor 
1/[1 — exp(27ria)], this factor being negligible when 9a is positive and large, as we assumed 
here. A similar conclusion holds when 9a, be negative and we replace a, z with ae +7 ™, ze +Vi . 



4. Numerical verification of the results. 

To verify numerically the uniform expansion we have used the real representation (3.11) for 
7* (—a, —z) with large positive values of a and z. We have used the recursion formula 

j* (-a - 1, -z) + zj*(-a, -z) = -- sin7rae z r(a + 1) (4.1) 

for testing the results. This relation easily follows from the well-known recursion aj(a, z) — 
j (a + l,z) = z a e~ z and the relation between 7(0, c) and j*(a,z) given in (1.2). To avoid 
overflow and strong oscillations we have normalized the function 7* (a, 2;) by using a quantity 
7 a (z) defined by 

7*(-a, -z) = z a cos7ra + sin7rar(a) e z j a (z). (4.2) 
The recursion reads for this new function reads 

-7„+i (z) + A % (z) + - = 0, A = -. (4.3) 
7T a 

Large quantities in the right-hand side of (3.11) are z a and z a exp(|a?y 2 ), and when dividing 
(3.11) by e z T(a) we can use the relation 

e z T(a) ~ V 2^ r*(a)' 
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where T*(a) is defined in (2.5). This gives 

la (?) 



irT*(a) 



I F {vyfi) + l a T a (rj) 



(4.4) 



Dawson's function F and T*(a) are computed with extended precision (about 19 relevant 
digits). We have used expansion (3.10) in (4.4) and truncated the series after the term C 6 (ry)/a 6 , 
with the expectation that the order of magnitude of the remainder is about 10~ 14 if a > 100. 

We have used the representations of the coefficients C ,...C 6 as given in Section 3 of 
Temme (1979). For < 1 we have used Maclaurin expansions of the coefficients. An 
extensive set of Maclaurin coefficients is given in Didonato &: Morris (1986). We have 
used enough Maclaurin coefficients in order to obtain 14 digits accuracy for Co on the interval 
\rj\ < 1, 12 digits for Ci, etc., with the intention to use the asymptotic expansion if a > 100. 

In the numerical verification we fixed a = 100 and took z in intervals around the point a, 
the transition point. In fact we used random numbers for z in intervals [(1 — 2 _fc )a, (1 + 2~ k )a], 
for k = 0, 1, ... 6. All computed left-hand sides of the recursion relation (4.3) were (in modulus) 
less than 5.5 x 10~ 18 . Values of 'ja(z) and 7*(— a, —z) around the turning point a = z = 100 
are given in Table 4.1, together with the error in the computation of ■~f a (z). 



z 


la (z) 


7*(-a, -z) 


error 


90.0 


2.464370784806670 


1.20552423411674426e+196 


5.3e- 


-18 


91.0 


2.372768700058633 


3.39249458130961502e+196 


2.7e- 


19 


92.0 


2.245263163760808 


9.45984032045937199e+196 


5.4e- 


20 


93.0 


2.081986332010084 


2.61247251775844538e+197 


7.0e- 


19 


94.0 


1.884270784695821 


7.14192476088520721e+197 


5.4e- 


20 


95.0 


1.654640974379185 


1.93183088802899244e+198 


3.8e- 


-19 


96.0 


1.396740233402602 


5.16774449595400219e+198 


1.4e- 


-18 


97.0 


1.115197442856446 


1.36641025407873601e+199 


4.9e- 


-19 


98.0 


0.815441855478639 


3.56901459711775166e+199 


l.le 


-19 


99.0 


0.503478153282224 


9.20223634675090428e+199 


1.2e 


-18 


100.0 


0.185636311520584 


2.34010604791689845e+200 


7.0e- 


19 


101.0 


-0.131687939165396 


5.86259361067072689e+200 


8.1e- 


19 


102.0 


-0.442286347630879 


1.44483553329971751e+201 


9.5e- 


19 


103.0 


-0.740370670716527 


3.49587466488947244e+201 


4.6e- 


19 


104.0 


-1.020772405557316 


8.28106485635440568e+201 


8.7e- 


19 


105.0 


-1.279100681803916 


1.91262956371154387e+202 


5.4e- 


20 


106.0 


-1.511852726452274 


4.28009100876763689e+202 


1.2e 


-18 


107.0 


-1.716474855286176 


9.18429865111133800e+202 


2.0e- 


-18 


108.0 


-1.891375273470379 


1.85454711082400020e+203 


4.9e- 


-19 


109.0 


-2.035892840388064 


3.38697191526660474e+203 


2.7e- 


-19 


110.0 


-2.150228198004453 


5.01390135464872193e+203 


l.Oe 


18 



Table 4.1. Values of ^ a {z) of (4.2) and 7*(— a, — z) for a = 100.25 with z— values around 
the transition point a = z; the "error" is the computed left-hand side of (4.3). Notation: 
1.20552423411674426e+196 means 1.20552423411674426 x 10 196 , 5.3e-18 means 5.3 x 10" 18 . 
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